Interaction of gravitational waves with matter 
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We develop a unified formalism for describing the interaction of gravitational waves with matter 
that clearly separates the effects of general relativity from those due to interactions in the matter. 
Using it, we derive a general expression for the dispersion of gravitational waves in matter in terms 
of correlation functions for the matter in flat spacetime. The self energy of a gravitational wave 
is shown to have contributions analogous to the paramagnetic and diamagnetic contributions to 
the self energy of an electromagnetic wave. We apply the formalism to some simple systems - free 
particles, an interacting scalar field, and a fermionic superfluid. 
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I. INTRODUCTION 

The interaction of gravitational waves with matter is 
important in a number of different contexts. One is in 
connection with the continuing quest to detect gravi- 
tational waves experimentally. Pioneering experiments 
were carried out using large metal bodies as detectors 
[l[ and this line of investigation has been further pur- 
sued at many centers. These detectors were designed to 
detect gravitational waves by observing the resonant ex- 
citation of elastic modes of the bars and the standard 
theory of such detectors uses elastic theory to calculate 
the response 

However, there are a number of suggestions that the 
response of matter could be very different from what is 
predicted on the basis of elasticity theory. One is that 
the direct coupling of the gravitational wave to electrons 
could enhance the absorption cross section @. For solids, 
the effect of electronic degrees of freedom has been taken 
into account within the Frohlich model for the electron- 
phonon interaction, which does not take into account ex- 
plicitly the long-range character of the Coulomb interac- 
tion, and the authors conclude that the effect of including 
electron degrees of freedom explicitly is very small Q . A 
more recent suggestion is that a superconducting metal 
would be a reflector of gravitational waves because ions 
and superconducting electrons respond in different ways 
to a gravitational wave, thereby creating a large electro- 
static energy that renders the superconductor "stiff" to 
the propagation of the wave These proposals under- 
score the need for a theory of the interaction of gravi- 
tational waves with matter that treats coupling of the 
gravitational wave to matter on a unified footing, takes 
into account the microscopic degrees of freedom, and also 
is able to include the effects of interactions. 

A second important area is the interaction of gravita- 
tional waves with astrophysical matter. Since much of 
this matter is diffuse and weakly interacting, the com- 
mon approach to this problem is to calculate trajectories 



of free particles in the presence of the curved spacetime 
produced by the gravitational wave. A review of early 
work on the dispersion of gravitational waves may be 
found in Ref. @. The effects of electromagnetic fields 
are included in some cases: for example, Servin, Brodin 
and Marklund [7, 8] showed that a magnetic field can ro- 
tate the polarization of a gravitational wave. Their work 
suggests that the role of the electromagnetic field is cru- 
cial to understanding the response of a charged system. 

The purpose of the present work is to develop a general 
formalism for describing the interaction of gravitational 
waves with matter. In particular, we wish to separate 
the effects of general relativity from those of calculat- 
ing correlations in the matter. Our approach is modeled 
on the semiclassical theory of interaction of electromag- 
netic fields with matter, in that we shall treat the grav- 
itational radiation classically. However, the matter will 
be treated quantum mechanically. The response of the 
system is calculated in a systematic way from a path- 
integral approach. We find that there are contributions 
to the response of matter to a gravitational wave that 
are analogous to the paramagnetic and diamagnetic re- 
sponses of a conductor to an electromagnetic field, and 
we give general expressions for them. Earlier work on 
interaction of gravitational waves with condensed matter 
d, Q has generally focused on the paramagnetic term, 
while in astrophysical applications the diamagnetic term 
often dominates. The formalism described in this article 
provides an economical way of deriving results for simple 
situations that have been considered earlier, while at the 
same time being of sufficient generality to be applicable 
to interacting many-body systems. 

The paper is organized as follows: in Sec|II]we develop 
the formalism for calculating the dispersion relation for 
a gravitational wave propagating in matter. Section IHII 
treats the case of free particles, both non-rclativistic and 
relativistic. In Sec. lIVI we analyze the coupling of a grav- 
itational wave to two interacting systems - a scalar boson 
field with a 0^ interaction (a Bose-Einstein condensate) 
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and a superfluid with paired fermions described by the 
Bardeen-Cooper-Schrieffer (BCS) theory. We describe 
possible directions for future research in Sec. |Vl 



includes the effects of both quantum-mechanical and sta- 
tistical averaging in the standard manner [T2j . The quan- 
tity Sm is the contribution to the action from matter, and 
it may be written in the form 



II. BASIC FORMALISM 

In a gravitational wave, the metric tensor 
gfi^{x.,t) deviates from the Minkowski metric 
Vfiiy — diag(l, —1, —1, —1) and we write 



5i^^(x, t) = ?7^^ 4- /i^,y(x, t) 



(1) 



where /ipi,(x, i) is the disturbance in the metric tensor. 
In keeping with the general approach we adopt, the met- 
ric tensor will be treated as a classical quantity. There 
is much freedom in the way in which the disturbance 
of the metric tensor is described and, for gravitational 
waves, a convenient choice is the transverse traceless 
(TT) gauge, in which h has only spatial components. 



d^h" = DM = 0, and K ^ 



For gravitational 



waves, the quantity /i^j^ plays a role similar to that of 
the potential in electromagnetic theory. In treating 
the effects of matter we shall assume that, in the absence 
of gravitational waves, space-time is flat. This is a good 
approximation for wave numbers small compared with 
the scale of the curvature tensor 0. 

To describe the interaction of gravitational waves with 
matter, we generalize to gravitational waves the semiclas- 
sical theory of electromagnetic response [ll| . The gravi- 
tational field is treated classically, but matter, including 
electromagnetic radiation, is treated quantum mechani- 
cally [13] ■ The purely gravitational contribution to the 
action is 
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dx„ dx'^ ' 



(2) 



where the second expression is the leading contribution 
for small h, k — SttG/c^ {G being the Newtonian gravi- 
tational constant), g = det R is the Ricci scalar and 
the integrals are taken over space and time, df^x. 

The contribution of matter to the effective action of the 
gravitational field due to matter is obtained by integrat- 
ing over all possible paths for the quantum mechanical 
motion and is given by 

5'cff(/i) =lnZ(/i). (3) 

Here the partition function Z{h) is given by [ll| 



Z{h) 



2?(V5,V',^)e"^'"^'^'^''^'''\ 



(4) 



where A is the electromagnetic potential and the fields 
tp and 'ip describe the other degrees of freedom of the 
matter. By taking the integration over the complex time 
coordinate to run from to i/3, where /3 = 1/T is the 
inverse temperature, one obtains a compact result which 



(5) 



where C is the Lagrangian function. The equation for 
the deviation of the metric tensor is found from the ex- 
tremum of the total effective action for the gravitational 
field, S'p-rav + '5'cff , and has the form 
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5lnZ 



(6) 



where □ 



is the d'Alembertian operator. 



This is equivalent to the standard result D/iy = —2KTi 
where Tij is the energy-momentum tensor, since 



_ SlnZ 



(7) 



In the TT gauge, h has no time components and therefore 
the indices run over the three spatial coordinates. For 
definiteness, we shall consider a plane gravitational wave 
propagating along the z axis, in which case the indices i 
and j can be either x or y. 

For /i — > 0, S In Z/Sh^^ is independent of time and is 
irrelevant so far as gravitational waves are concerned. 
Expanding SlnZ /Shij to first order in h one finds 



where 
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(9) 



Thus the quantity T,fj plays the role of a self energy 
for the gravitational wave. For simplicity, we shall con- 
sider a medium that is isotropic, spatially homogeneous, 
and invariant under time reversal. It is then convenient 
to work with the quantities h+ = {hxx — hyy)/2 and 
hx = hxy = hyx that correspond to normal modes of the 
system, and these satisfy the equations 



and D/ia 



0, (10) 



the factor of two being due to the fact that in Eq. ([5]) kl 
can be both xy and yx. 

In order to find an expression for S for small pertur- 
bation in the metric, we expand the Lagrangian density 
in powers of h, 




h=0 



2 Sh'^Shki 



h'^ hki + Oih^) . (11) 
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Quite generally, the energy-momentum or stress tensor 
is given by [1^ 



and therefore 



h=0 



-g 5W 



2 



(12) 



(13) 



where Tj^- is the stress tensor in flat spacetime. 

Since the self energy is a second functional deriva- 
tive with respect to h, it is important that all quanti- 
ties of second order in h are calculated consistently. In 
particular, care must be taken to distinguish upper and 
lower indices, as one can see by noting that the condition 
ff*'^ gkj = implies that 



(14) 



to second order in h. Moreover, remembering that gij 
^tj f^ik o{h'^), one finds 



Ttj - gucT gij 



(15) 



After some calculation, the equation for the field is ob- 
tained: 



ah, 
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where 




(16) 



(17) 



is what we shall refer to as the "diamagnetic" contribu- 
tion and, with arguments written out explicitly, 

mr^^)) dv^dt^9{t,^t^){[T^^{l),TU2)])hu(2) 

(18) 

is the "paramagnetic" contribution. Here (O) = 
/ D[4j, V', A] Oe-'5'"('^''^'^'''=o) denotes the thermal aver- 
age of the operator O in flat spacetime, [A, B] denotes 
the commutator, and 9{t) is the unit step function. If 
/ifci varies in time as e~"^*, one therefore finds 



(^TP--(ri,t)> 



where 

X^J,/c^(rl,r2) 



c'r2Xu,fci(ri,r2)(5/ifei(r2), (19) 
''(i;"(ri))„™(T0(r2))™„ 



En - E„i + Lo + irj 
(T0(r2))„™(T0(ri))„„l 



En — E„ 



LO — 17/ 



(20) 



where P„ is the statistical weight of energy eigenstate n 
for the system in fiat space. 



If the Lagrangian density is local in time, the diamag- 
netic term is independent of frequency and therefore be- 
haves as a mass term. For matter described using a non- 
relativistic framework, the Lagrangian density may be 
nonlocal in space, in which case the mass will depend 
on the wavevector. The paramagnetic term depends on 
both frequency and wave vector and contains information 
about excited states of the matter. 

We now examine the various contributions to the re- 
sponse of the stress tensor. The first term on the right- 
hand side of Eq. (|T7)) is proportional to the pressure of the 
matter and it is therefore independent of the frequency 
of the gravitational wave. If the Lagrangian density is 
local in time (as is usually the case), the second term is 
also frequency-independent. Thus both these terms be- 
have as a mass term for the gravitational wave, which is 
why we refer to their total as the diamagnetic term. 

Notice that the expectation value of the second deriva- 
tive of the Lagrangian in vanishes identically for free 
particles, leaving only terms proportional to the pres- 
sure; this is in accord with the previous results of the 
relativistic literature [1, |^. However, this term gives a 
contribution to the dispersion relation in the presence of 
an interparticle potential. For example, this contribution 
must be taken in account for calculating the self energy 
of a gravitational wave interacting with a scalar field (a 
Bose-Einstein condensate), as we shown in Sec. IIV Al 

The only contribution to the self energy that has 
an energy denominator and can thereby give rise to 
absorption of gravitational waves by matter is the 
paramagnetic term. The expression ([T5| is analogous to 
that for an electromagnetic wave in a conductor, where 
the response of the matter is evaluated self-consistently 
using linear response theory. The paramagnetic term 
does not contain the current-current response, but an 
equivalent expression with the stress tensor. For an elec- 
tromagnetic wave the diamagnetic term is proportional 
to the particle density, while for a gravitational wave it 
contains terms proportional to the pressure and terms 
proportional to the second derivative of the Lagrangian 
density. In the equation for the field hij the right hand 
side contains only quantum mechanical and thermal 
averages of quantities in fiat space. Thus in the present 
approach, the effects of general relativity have been 
decoupled from the problem of solving the many-body 
problem for the matter. This is possible because the 
gravitational fields are weak. 

A number of properties of response functions at long- 
wavelength may be obtained by a consideration of con- 
servation laws. This has previously been done for the 
density, spin density and current responses in the con- 
text of Fermi liquid theory [343 and we here apply 
these ideas to the stress-tensor response. 

We consider the response of an initially uniform 
medium to the application of a perturbation having the 
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form 



(21) 



where J7_q is the strength of an applied external field and 
"h. c." denotes the Hermitian conjugate. If the operator 
O satisfies a local conservation law, in coordinate space 
the operator equation for the conservation law is 



(22) 



where jo is the operator for the corresponding current. 
On taking matrix elements of the Fourier transform of 
this relation between energy eigenstates of the unper- 
turbed system, which are labelled by m and n, one finds 



^(a 



q-(j?)n 



(23) 



This equation demonstrates that, if O satisfies a local 
conservation law, and provided the corresponding cur- 
rent is not divergent in the long- wavelength limit, then 
matrix elements of (C'q)^™ vanish for q — )■ for all states 
for which cj„,„ ^ 0. Expressed in other terms, this states 
that the only nonvanishing matrix elements of (C'q)„m 
are between states whose energy difference falls off at 
least as rapidly as q. It is this observation when applied 
to the particle density, and in the case of translationally 
invariant systems also the particle current density, that 
lies behind the success of Landau Fermi liquid theory in 
providing a powerful way of parametrizing the proper- 
ties of long-wavelength properties of normal Fermi liq- 
uids. The calculations of response functions for interact- 
ing systems in Sec. IIVI will illustrate these general prop- 
erties, but first we describe results for free particles. In 
this case there are matrix elements of the stress tensor 
operator to states having nonzero excitation energy van- 
ish in the long-wavelength limit because the momentum 
of a particle and its velocity are both conserved quanti- 
ties and, consequently, the contribution of a particle to 
the stress tensor is conserved. 



III. FREE PARTICLES 

In this section, we study the response of a system of 
noninteracting particles to a gravitational wave. First, 
we consider free particles and derive simply results pre- 
viously obtained by other methods in the astrophysical 
literature. In addition, we explore two other systems 
where the quantum mechanical nature of the system is 
relevant: the Bogoliubov theory of a Bose-Einstein con- 
densate and the ECS theory of superfluid fermions. 

We begin by considering the case of a noninteracting, 
nonrelativistic particles. The diamagnetic contribution 
to the self energy of the gravitational wave, Eq. pT|) 
may be calculated simply for free particles obeying either 
the Schrodinger equation or the Klein-Gordon equation, 
since (£) is zero and therefore {S^ ^ — g C / Sh"^) = in 



the TT gauge, and the diamagnetic contribution con- 
tains only the pressure P = (Ta); as we shall show, for 
nonrelativistic particles the diamagnetic contribution is 
smaller by a factor (w^)/c^, where (u^) is the mean square 
particle velocity. Thus one finds 



2 „2 



327rGP 



c q 



(24) 



to first order G and first order in (w^)/c^. For fermions 
at zero temperature the pressure is P = npjp/5m, where 
Pi? is the Fermi momentum; the dispersion relation is 
therefore 



„2 „2 



c q + 



327r 



Gr 



(25) 



For an ideal gas obeying Maxwell-Boltzmann statistics, 
the pressure is P = nT, and the dispersion relation be- 
comes 



o;^ w q^ + 32n Gmn- 



T 



(26) 



We see that the dispersion relation depends on the 
"Jeans" frequency = {GnmY/'^ characteristic of 
gravitational collapse, reduced by a factor (w^)^/^/c. The 
result (1^^ coincides with the result in the literature 

BE 20]. 

We turn now to the paramagnetic contribution. In or- 
der to simplify the discussion, it is convenient to consider 
the specific case of a gravitational wave corresponding to 
a disturbance of h^y = hyx = hx propagating in the z- 
direction. The stress tensor for free particles in flat space 
is given by 

T%(l) = -^(Vi-ViO.(Vi-ViO, ^t(i')^(i) 

4 m I'^i 

(27) 

where "0^(1) is the particle creation operator and "0(1) the 
annihilation operator at the point (ri,ii). In a uniform 
medium, it is convenient to work with the spatial Fourier 
transform of this quantity, which is given by 



(TO 



qjxy 



E PxPy 
«p-q/2«P+q/2 

P 



(28) 



In the literature the response to a gravitational wave 
is often studied by using the Vlasov equation to model 
the behavior of the excitations in the system. In this 
section we first present the Vlasov equation approach and 
then show that the same results may be obtained by a 
quantum mechanical treatment. 

In a homogeneous system, the perturbation in the 
stress tensor of a non-interacting gas is given by 

5TP-'"(q) = |^t^^(p)^np(q), (29) 



where 



(30) 
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is the stress tensor associated with a single particle, and 
the Vlasov equation reads 



d 



Snp + v VrSnp - (Vrt° h^y) ■ Vpn° = . (31) 



dt 

Thus one finds 
where 

Xoi:y,xy — ^ 



dp pIpI q-Vpn° 



{2t:Y fTT-^ w — q-p/m 



(32) 



(33) 



is the transverse stress-tensor-stress-tensor response 
function. 

For gravitational waves, the frequencies of interest are 
approximately cq and therefore for nonrelativistic parti- 
cles one may expand the denominator in Eq. p3p and to 
leading order in the result is 



JTPr'' ^ ^ 



1 f dp pIpI (q-p)2 9n° 



xy 



ij^ J (2n)^ m? m? dtp 



•'xy ■ 



(34) 



For nondegenerate particles, the distribution function is 
Maxwellian and one finds 



Xxy,xy - " ^ ^2 



while for a Fermi gas at zero temperature 



1 



Xxy,x'y 
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35 ^a;2 



(35) 



(36) 



where in Eq. ([36)1 we used the fact that n = N/V = 
lyp^/dir^, where u is the number of degenerate inter- 
nal states of the particle, due to spin, isospin or other 
symmetries. By including the paramagnetic term, the 
dispersion relation becomes 



2 2 I 

c q + 
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167rG 



Xxy,xy 



(37) 



The results may also be obtained from a quantum- 
mechanical calculation based on Eq. (15^ . In the nota- 
tion of second quantization, the stress tensor operator for 
a free particle system is given by 



f,y{i) = (Vi-VlOx (Vi-ViO^ V'^(i') ^(1) 

(38) 

where to simplify the notation we do not write sums over 
internal states explicitly. The expectation value of the 
stress tensor therefore reads 



(f,,(i)) = -— (Vi-ViO.(Vi-ViO,e(i,i')li'^i+ , 

(39) 

where C/(l,l') = — {il}'^ {!') is the (finite tempera- 

ture) single-particle Green function. The paramagnetic 
response is given by 



(5f P""(l)) - i / dr2dt2 (Vi - ViOx (Vi - ViOy (V2 - VaOx (V2 - V 
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■2')y 



V;t(l')V^(l),^^(2')^(2)]) h,y{2)., 

2' ^ 2+ 



which in Fourier space becomes 

(5fPp(q,ic^„)) = -Y,^-^QiP + Ci^UJl+lUJn)Q{p^UJl)hxy{ci^LOn) 



9 *? 
PxPy 



Plpln°{p + q)-n^{p) , , . N 



m? i cj„ - Ep+q + Ep 



(40) 



(41) 



where the Matsubara frequencies are a;„ — 2'KnT and 
LOi — 21ttT for bosons {{21 + 1)ttT for fermions). For 
q <^ (p) , Eq. (|1T|) reduces to Eq. ([55)) obtained from the 
Vlasov equation. 

We now comment briefly on the case of relativistic 



particles, and for definiteness we shall consider particles 
described by the Klein-Gordon equation. As remarked 
above, for such particles the diamagnetic response is 
given in terms of the pressure, just as for nonrelativistic 
particles. The calculations for the paramagnetic response 
may be performed essentially as before and for q <^ (p) 
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the effect is to replace the mass m by the "relativistic 
mass" m(l + (p/mc)^)-^/^. In general, one cannot assume 
that the particle velocity is small compared with c, and 
consequently the response function has to be evaluated 
numerically. Simple results may be obtained for ultrarel- 
ativistic particles, since the particle velocity is c for all 
momenta. In that case, the integrals over the polar angle 
and the magnitude of the momentum decouple and one 
finds 



Xoi:y,xy 



dp 



-P( 



16 
15 



2s'' + s{s 



s — II 
10 2 

2 



(27r)3 ~dp 



1)2 In 



s - 1 



(42) 



where p — p/p, p — p ■ 2, and s = uj/cq. The response 
function for the transverse components of the stress ten- 
sor is finite for s = 1, unlike the density and current 
response functions, which have a logarithmic divergence. 
The physical reason for this is that particles moving in 
the direction of propagation of the wave give vanishing 
contributions to the transverse components of the stress 
tensor. For s — 1, Xxy,xy = 8P/45, which is of the same 
order as the diamagnetic response. This result can be 
added to the diamagnetic contribution, Eq. ([M)) . to give 
the following dispersion relation for the ultrarelativistic 
case: 



2 2 

c q 



45c^ 



P. 



IV. INTERACTING SYSTEMS 



where, in the Gross-Pitaevskii approach, ?7o = Airfi^a/m, 
a being the scattering length for two-body scattering, is 
the strength of the effective two-body interaction. There- 
fore the spatial components of the stress tensor (fT^ are 



1 



T,,(x) = -_(Vi-ViO.(Vi-ViO, ^^lOV^i) 



(45) 



The last term in ps)) is the pressure due to the in- 
terparticle interaction and in the Bogoliubov approxima- 
tion, in which is replaced by a c-number ^/n^ with hq 
the condensate density, it becomes 71^1/0/2. If depletion 
of the condensate may be neglected, no may be replaced 
by n and the result agrees with the one obtained from 
the thermodynamic relation P = n^diE /n)/dn, where 
E is the energy density. 

The interaction does not contribute to the paramag- 
netic term in the dispersion relation, because the grav- 
itational wave is transverse, but there is a diamagnetic 
term since 




(46) 



where P ^ Uq 12^/2 is the pressure. The contribution 
^ydia gjvgg then the dispersion relation 

nUo 



2 2 
c q 



IG-kGi 



(47) 



We now consider the paramagnetic term, which is not 
generally zero. In Fourier space the contribution to the 
stress tensor from the kinetic energy may be written as 



q. . 

-)j flp+q 



In this section we consider two examples of interacting 
systems at zero temperature. The first is an interacting 
boson field. Due to the interaction, this has a nonvanish- 
ing diamagnetic contribution to the response to a gravita- 
tional wave. Both it and the BCS superfiuid have param- 
agnetic contributions due to excitation of two excitations 
with nonzero energy even at long wavelengths, and serve 
as an illustration of the general results described at the 
end of Sec. Ill 



A. Interacting boson field 

The Lagrange function for a nonrelativistic boson field 
(a Bose-Einstein condensate) with a short-range interac- 
tion is 

2m 2 




(44) 



-. — Vrkiqi qj (flq + a_c 
4to 



1 

TO 



P 

Pi Pj al flp+q , 

p 



(48) 



where the prime in the sum over p denotes that only 
the atoms outside the condensate are taken in consider- 
ation. The condensate contribution vanishes identically, 
because the only non-zero component of q is in the z di- 
rection. This component never appears in the transverse 
response, i.e. qi = qj = for i,j = 1, 2. 

Elementary excitations of the condensate are created 
by operators cij^ and destroyed by cik, which are related 
to the the particle creation and annihilation operators by 
the Bogoliubov transformation 



flk = Wk - Wk alk , 
dl ^uual- Wk a-k 



(49) 
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with ul = l + vl = [l + {eu + n (7o)/a;k]/2, ek = m 
and = + 2nUo Ck- The stress tensor is given by Eq. 
psp and the only contribution that gives nonzero matrix 
elements when acting on the ground state is that which 
creates two excitations, which is given by 

(50) 



where we have also made use of the fact that the gravi- 
tational wave is transverse, and therefore — % = 0. 

Inserting the expressions above for the matrix elements 
of the stress tensor operator and the excitation energies 
into the general result pOI) . one finds at zero tempera- 
ture, and when analytically continued to a real frequency 
w, the result 



dp 

(2^ 



2 9 
PxPy , 



l^p+q Up) 



+ 17]- [Ujp+q +UJp) + IV + i^p+q + ^p) - 



(51) 



Equation ((5T|) shows that the response does not vanish 
even for q = 0. It is not possible for a gravitational 
wave, which is transverse, to excite a single Bogoliubov 
excitation because the latter is longitudinal. It is rel- 
evant to stress that, while the contribution from states 
with two excitations to the transverse current-current re- 
sponse function vanishes in the long wavelength limit, the 
corresponding contribution to the stress-tensor-stress- 
tensor response remains nonzero, due to a different sign 
in the Bogoliubov factors inside the parenthesis in ([51]) . 
This is a specific example of the general result given in 
Sec. In] and is a consequence of the fact that the stress 
tensor does not obey a conservation law when there are 
interactions. 

The integral in Eq. (ISTj) is ultraviolet divergent. This is 
due to the fact that we have used an effective low energy 
theory to calculate a quantity that cannot be expressed 
in terms of the constants in the theory. However, the 
response at frequencies with a magnitude of order nUo or 
less may found by subtracting from the response function 
its zero frequency value. For simplicity we consider the 
long- wavelength limit, g — ^ and find 



Xxy,xyi<l,i^) - Xxy,xy{ci,0) ^ 



(2up Up) 



dp pIpI 



(27r)3 m2 



2Wr 



2wc 



(52) 



The integral may be performed analytically, but the re- 
sult gives little insight. The imaginary part is simple, 
and is given by 

1 {nU^Y p{u:f 
Inix(q,w) = ^ 7, rr , 9/ (53) 



where 



1577^ ui 2nmUo + (lu) ^ 

1/2 



p{uj) = 



(54) 



At frequencies much less than nU^ this varies as w** while 
for frequencies much larger than nJ7o (but still small 
enough for the low-energy theory to be valid) it varies 



1/2 



as w 



B. Fermionic superfluid 

In this section we compute the response, at zero tem- 
perature, of a superfluid made up of pairs of fermions in 
two internal states, which we refer to as up and down. 
The contribution to the stress tensor operator from the 
kinetic energy may be written as 



X (ap+qt«Pt + «-(p+q)ia-Pi)-(55) 

There will generally be in addition a contribution from 
the interaction energy but we ignore this since, for weak 
coupling, it is small compared with that from the kinetic 
energy. We shall assume the superfluid to be of the BCS 
type, with s-wave pairing between two spin states, de- 
noted by "up" and "down" . The elementary excitations 
are linear combinations of particles and holes, which are 
destroyed by operators dk, /3k and created by the Her- 
mitian conjugate operators. The particle creation and 
annihilation operators are related to these by the canon- 
ical transformation 



Wk,5-1. 



k ' 



(56) 



where u\ = \-v\ = 1/2 (1 + ^q/Eq), £.q = eq - f-i, and 

Eq = y^A^ + Q is the energy of an excitation, with A 

being the energy gap. For simplicity, we restrict ourselves 
to zero temperature, and therefore the contribution to 
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the stress tensor operator coming from the kinetic energy 
is 

TtjicO = ^ + |)» {p + |)j {Up+q Vp + Vp+q Up) 

P 

-(p+q)"p)-' (57) 

where we have neglected terms like $ and a (3^ that 
vanish when acting on states with no excitations present. 
For weak coupling, the contributions to the stress tensor 
operator from the interaction energy will be small com- 
pared with those from the kinetic energy, and we shall 
neglect them. 

As before, the diamagnetic term is proportional to the 
pressure of the fluid. In this section we consider the first 
term of the dispersion to be the same as the one for the 
free fermion case, by neglecting the effect of the interac- 
tion on the pressure. Therefore, the dispersion relation 
reads 

2 2 2 327r^ v% 167rG 

uj M c q + ——Gmn— H — Xxy,xy ■ (58) 

o c c 

The response function can be found by inserting the ex- 
pressions for the matrix elements into the general expres- 
sion (Uni) 

Xxy,xy[ti,^) = 2^ ^ [Up+qVp+Vp+qUp) 

P 



Lo - {Ep+q + Ep) +iri 



UJ + {Ep+q + Ep) +17]^ 

m2 ^ ^''P' Ep{{co + iv)^-AE^)' 

(59) 



where the final expression is valid for q ^ A/vp. The 
imaginary term is nonzero only for a; > 2 A, since the 
minimum energy of a single excitation is A. 



The sum in (j59p gives a large contribution for ^ 2 A, 
and far from this resonant frequency the response hardly 
makes any contribution, since the self energy is propor- 
tional to G/c^. More quantitatively, this sum can be 
transformed to an integral, and this integral can be solved 
in spherical coordinates. Integrating over the angles we 
have 



X.....(q,-) = J^J^^l dpp 

A2 

Ep{{u; + t^)^-4E'^) • 

(60) 



This integral does not converge, the problem being that 
an effective low-energy theory has been used to calculate 
a quantity that has important contributions from high- 
energy states. To investigate the low-frequency structure 
of the response function, we subtract from the response 
function, its value for lo ~ 0. The integral converges and 
consequently one may evaluate it putting A equal to its 
value at the Fermi surface. 



Xxv,xv{(l,i^) - Xxv,xv{(i,UJ = 0) = -fr „ 2^ 2 / ^PP^ 



i?3((.. + i^)2_4£;2) 
1 f°° 1 



5 J^oo ^(^2+A2)3/2((c^ + i^)2_4(^2 + A2)) 

= lnpFVFF{uj/2A) , (61) 
5 

I 

where, for a; < 2 A with ft = sine/), and for a; > 2 A 

F{n) = i + v/f^2~~T) ^ 



sm2(h ' The real and imaginary part of F are plotted in Fig. [T] 
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FIG. 1: Behavior of the function F{uj,A), Eqs. (I62| and (|63ll . 
The full line is the real part of F and the dashed line is the 
imaginary part. 



This calculations represents the simplest approxima- 
tion for the response, but they do not take into account 
residual interactions between excitations. Such interac- 
tions are important for the collective behavior and lead, 
e.g., to the Bogoliubov-Anderson sound mode [13, HI], 
which represents a density wave in the condensate. How- 
ever, the effects of the residual interaction on the stress- 
tensor response function are expected to be less dramatic 
since the it is transverse, not longitudinal. The present 
case is more analogous to excitons in superconductors 
24 1 and pairing vibrations in atomic nuclei (251, where 
the momentum dependence of the interaction plays a cru- 
cial role. 



V. CONCLUDING REMARKS 

In this paper, we have developed a general framework 
for studying the interaction of a weak gravitational wave 



with matter. A virtue of the approach is that it sep- 
arately clearly the effects of general relativity from the 
problem of solving the many-body problem for the mat- 
ter. The matter gives a self energy to the propagator of 
the gravitational wave. This self energy has contributions 
analogous to the paramagnetic and diamagnetic contri- 
butions to the self energy of an electromagnetic wave in 
matter. The contribution corresponding to the param- 
agnetic term is proportional to the stress-tensor-stress- 
tensor correlation function for the matter. Because the 
stress tensor operator is not a conserved quantity, ex- 
cept for non-interacting particles, this correlation func- 
tion does not in general vanish in the long wavelength 
limit for nonzero frequency, and we illustrated this by 
explicit calculations for a Bose-Einstein condensate and 
a BCS superfluid. The general formalism in this paper 
makes for a very simple derivation of the dispersion rela- 
tion for gravitational waves in astrophysical plasmas. 

There are a number of possible directions for future 
work. In this paper we have considered only an infinite 
medium, and one could extend the treatment to take into 
account the effect of boundaries. Another application is 
to systems, like metals and superconductors, in which 
the Coulomb interaction plays a key role. The formalism 
may also be used to establish the relationship between, 
on the one hand, the microscopic theory in terms of par- 
ticles and their interactions and, on the other hand, elas- 
tic theory which has been commonly used to discuss the 
response of gravitational wave antennas. 
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